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Abstract 

We investigate a simplified model of dark matter where a Majorana fermion x coannihilates 
with a colored scalar top partner t. We explore the cosmological history, with particular emphasis 
on the most relevant low-energy parameters: the mass splitting between the dark matter and the 
coannihilator, and the Yukawa coupling y x that connects these fields to the Standard Model top 
quarks. We also allow a free quartic coupling A/, between a pair of Higgs bosons and t pairs. We 
pay special attention to the case where the values take on those expected where t corresponds to 
the superpartner of the right-handed top, and x is a bino. Direct detection, indirect detection, and 
colliders are complementary probes of this simple model. 
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I. INTRODUCTION 


The weakly interacting massive particle (WIMP) is a theoretically attractive candidate 
for the dark matter, but direct detection bounds are becoming increasingly strong. Indeed, 
results from the LUX experiment [T] exclude, depending on the dark matter mass, spin 
independent dark matter-nucleon scattering cross sections above 10 -45 cm 2 . Such limits 
represent a challenge to dark matter whose scattering strength is related to annihilation 
via simple crossing symmetries. That is, tension arises if the dark matter thermal relic 
abundance is set by the diagrams that control scattering For example, this places 

pressure on so-called “well-tempered” scenarios [1] of the minimal supersymmetric Stan¬ 
dard Model (MSSM), wherein the proper dark matter abundance is achieved by titrating 
a non-interacting bino with the appropriate amount of Higgsino, which has full-strength 
annihilations under the weak interactions. 

There are several well-known examples where the crossing symmetry is spoiled [5]. One 
possibility is that a state Y co-inhabits the thermal bath with the dark matter x at the time 
of freeze-out, denoted by the temperature T F . This is possible if the mass splitting between 
X and Y is not too large, A M x y Tp. In this case, the dark matter abundance is determined 
not only by the size of the annihilation cross section of x 1° Standard Model (SM) fields, 
a(xx SM), but also by cr(xY —>■ SM) and, as long as conversions between x and Y 
are sufficiently rapid, a(YY —> SM). In the last two coannihilation cases, the crossing 
symmetry with direct detection is clearly spoiled: the interactions of the Y particles with 
the SM, vitally important for relic density considerations, are irrelevant for direct detection 
(at tree level). 

Given the current absence of direct detection signals, we are motivated to further consider 
the coannihilation scenario, with an eye towards elucidating the physics that determines the 
relic density. We consider a simplified model that captures important features of the case in 
the MSSM where a single light stop squark (say, tp) coannihilates with a bino-like neutralino. 

Colored particles make for interesting coannihilators since their large QCD cross section 
somewhat mitigates the otherwise strong mass degeneracy needed for a substantial correction 
to the relic abundance. Furthermore, the LHC clearly calls for a detailed consideration of 
possible scenarios which include strongly interacting new physics. For the particular case 
where the coannihilator is a top partner, direct detection constraints are less severe than 
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if the coannihilator has large Yukawa coupling with the first two generations, as might be 
expected for other “quark partners” [Hj. In additions, light top partners could provide a 
testable mechanism for electroweak baryogenesis 13 i- Finally, we note that it might be 
expected that top partners may be near the bottom of the spectrum of any new TeV scale 
physics if they are relevant to cutting off the quadratic divergence in the Higgs boson mass 
squared parameter. 

Much of the relevant computations of the underlying physics for stop coannihilation in 
the full MSSM were already done in the pioneering work of Ref. [2] and were recently 
reevaluated in HD]. Often, the MSSM phenomenology is discussed within the context of 
the Constrained Minimal Supersymmetric Standard Model (CMSSM), and it is not always 
transparent how to translate results from the m o - m 1/2 plane to the processes that underlie 
the determination of the relic density. And yet, there are often a very small number of 
processes that capture most of the early universe cosmology which can be understood in 
terms of a few low-energy parameters. This motivates a simplified model approach to the 
dark matter. Only the lightest states need be considered for computations of cosmological 
history, direct detection, and collider signatures. Steps in this direction for a theory of a 
dark matter accompanied by a colored partner were taken in PUMH]. If desired, a given 
UV model can be mapped on to this simplified model. 

This article is organized as follows. In Sec. |TT] we define our model and terminology. 


We then analyze the relic density in Sec. Ill This includes a basic review of the thermal 
evolution of a dark matter particle, including coaanihilations. Relevant cross-sections are 
presented with particular emphasis on their scalings with the parameters of our model. We 
also detail the effects of the Sommerfeld enhancement. The resulting consistent parameter 


space is then mapped out, specializing also to the supersymmeteric case. Sec. [TV] discusses 
current and future experimental probes of our model, including direct and indirect dark 
matter detection as well as collider experiments. We reserve Sec. [V] for our conclusions. 


II. EFFECTIVE STOP 

We consider a generalization of the neutralino stop coannihilation scenario. Our “stop” is 
a generic colored scalar with an arbitrary coupling y x to the top quark and the dark matter, 
which we take to be a gauge singlet Majorana fermion y. We consider the scalar in the 
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fundamental representation of the SU(3)c to be a partner of the right-handed top, so 


£'x = VxtRXt + h. c. , 


( 1 ) 


corresponds to the Lagrangian involving the dark matter held x interactions with SM 
fermions. The MSSM, with x identified with the bino, and t identified with the right- 


| \f2g Y -, with () Y the hypercharge gauge 


,MSSM 

X 


handed stop, corresponds to the case y: 


coupling. We also take the quartic coupling between the SM Higgs boson and a pair of stops 
A h to be free. The interaction of t with the SM is given by 


4 nt = \DJ\ 2 + \ h h ] htH. 


( 2 ) 


In the MSSM without stop mixing, \^ ISSM — y^ — 2 ntf/v 2 ~ 1, where v = 246 GeV and 
y t is the SM Yukawa coupling (we have taken the decoupling limit for the Higgs boson and 
neglected corrections from supersymmetry breaking). 

In the MSSM, allowing for stop-mixing with large H-terms - as might be motivated by 
the observed large value of the Higgs boson mass - can lead to interesting cases with large 
couplings of (both) stop states to the Higgs boson. The cosmology is potentially modified 
in interesting ways with respect to the simpler case presented here, and we plan to explore 
this in detail in an up-coming publication ra. 

It is possible that the scenario defined via Eqs. (JTJ) and (J2| could be realized in a simple 
extension of the MSSM by enlarging the gauge sector and identifying the LSP primarily 
with the new gaugino. If the corresponding H-terms do not decouple, it is also possible to 
modify the Higgs boson mass. However, our interest is not in model-building, but rather 
we will use this set-up as an effective parameterization of a model with a small number of 
degrees of freedom. 

III. RELIC DENSITY 

The observed relic abundance, f 1h 2 = 0.11805 ± 0.0031 SB. can be easily achieved by 
thermal freeze-out of weak-scale dark matter, an observation often dubbed the “WIMP 
miracle.” The dark matter abundance can be calculated to good accuracy as 



( 3 ) 


4 



with g* ~ 80 is the number of relativistic degrees of freedom and 

J(xf)= [ ^Y~dx, (4) 

Jx F x 

where a is the annihilation cross-section, v is the relative velocity and the thermally averaged 
cross section is defined as 

x 3/2 9 , 

W) = 2V75 1 (ar)v V dV ■ (5) 

The freeze-out occurring at x F = m x /Tp is determined by the iterative equation 

Xp = log 

where g — 2 for a Majorana fermion. In the absence of coannihilations, the thermally 
averaged annihilation cross section can be expanded as: 

(^xx v ) = MXX ->• SM)v) = a + b(v 2 ) + 0((v 4 )) = a + • ( 7 ) 


4.64 x 10 17 GeV g m x (<7v) 


1/21/2 


yielding 


Q x h 2 ~ 


5.77 x 10- n GeV- 2 a; F 

\J9*(xp ) T 36/iCiT’) 


As a rough rule of thumb the freezeout temperature is Tp ~ m x / 25, corresponding to an 
op? « 25; therefore 

n ^-012^ ftr’h a + 3fc/x F y 1 

» UrJ' soJ \3 x 10 -26 cm 3 /s/ ' ' 

In the model we consider, the dark matter particles would predominantly annihilate into 
a pair of tops, (xx ~t tt)i via the t and w-channel exchange of a t. The .s-wave contribution 
to this is given by: 

3 m 2 y x \ vn 2 — m 2 

a= W 3 - J~2 - 2, 2 1 2 ‘ ( 10 ) 

32 7rm x (m| — mf + m x ) / 

As a chirality flip of the top in the final state is required, a is proportional to m 2 , and 
therefore in the limit where x i s much heavier than the top, the dominant contribution is 
due to the velocity suppressed p -wave contribution: 


b „ m lVx ( m i + m x) 
167T (m 2 + m£) 4 
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We present only the rrit —>• 0 limit of b here, however the full mass dependence is always used 
in our numerical results, which are obtained by solving the relevant Boltzman equation(s) 
numerically with the help of MicrOmegas 3.3 [22], 

As we will review below, the above expressions must be modified if additional degrees 
of freedom, e.g. t, are still present in the thermal bath during freeze-out. In addition, it is 
well known that the annihilation cross section of charged non-relativistic particles can be 
modified by non-perturbative corrections: the Sommerfeld effect [23H26] . As t is charged 
under QCD, these corrections can become important and lead to significant modifications 
of the relic density Finally, we caution the reader that higher order corrections 

to stop coannihilations can be significant [m EH, and the accuracy of the relic density 
calculation lags behind the precision of the observations by Planck. 

A. Coannihilation 

The canonical calculation of the relic density needs to be modified if other particles of the 
dark sector are close in mass with the dark matter [5]. Not only dark matter annihilations, 
but also processes involving the next to lightest particles in the dark sector can contribute to 
an effective annihilation rate. Quantitatively, this corresponds to replacing the annihilation 
cross section er xx v in the Boltzmann equation with an effective cross section given by [5] 

eq eq 

n i n j / x 

W//V = 2^ n °*) 2 ^i v ’ ( 12 ) 

2 j ' < ^ k k ) 

where n\ = ^(m i T/(27r)) 3//2 e _mi / T ; vtii is the mass of the particle i, and r/, counts the 
number of internal degrees of freedom. The relic abundance can be approximately found 
using Eqs. now utilizing this effective cross section. We emphasize, however, that 

our numerical results rely on numerical solution to the Boltzman equation via MicrOmegas 
3.3 [22]. 

In the model studied here, coannihilation of x with t as well as the annihilation of 
tt* pairs are relevant for the calculation of cr e ffV pj. In general, the contribution due to 
processes with t in the initial state require a convolution of the annihilation cross section 
and its thermal abundance. The abundance of the heavier state is suppressed relative to 
that of y by additional factors of e _Am / TF , where Am = mi — m x . Due to this exponential 
dependence, the relic density is extremely sensitive to the mass splitting between m x and 
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rrif. Contributions to the relic density from the coannihilations of a pair of ts are doubly 
exponentially suppressed compared to yy annihilations, and the annihilations of x with 
a t are singly exponentially suppressed. For mi > 1.2 m x , coannihilations can safely be 
neglected. Note that the possibility of coannihilations provides a lower limit on Am as a 
function of m x . If the mass splitting is too small, then coannihilations will be too effective, 
owing to the large irreducible (QCD) cross section for tt* annihilations: 


7 q 4 

av(tt* ->• gg) = 


216 7rm? 


(13) 


In addition to {tt* —y gg), channels involving electroweak bosons can be relevant, spe¬ 
cially in the large A h regime. Because a relatively small number of processes contribute 
to the determination of the relic density in our model, we find it instructive to reproduce 
simple expressions for the relevant cross sections, tabulated in Table [T} Only the dominant 
contributions in the limit m x and rrif much larger than the weak scale, are listed. We have 
separated out the cross-sections in Table [T] into three parts, depending on the initial state. 
The relevance of each initial state depends on the mass splitting between x and t , which we 
parameterize by the ratio of their masses, r = mi/m x . All cross-sections were computed by 
implementing our Simplified Model in CalcHEP 3.4 [32]. 

As seen in Table [T] the channels scale differently with the dark sector couplings y x and 
A h- Cross sections of the processes initiated by y£, which are thermally less suppressed than 
tt *, are proportional to y x and therefore depend rather strongly on the dynamics of the new 
sector. Typically, the dominant channel here is (yt —> gt ), which scales as y x g 2 s , but, final 
states with the Higgs or electroweak gauge bosons are not irrelevant and the ht final state 
can receive a significant enhancement if A h > 1. On the other hand, {tt* —» gg), typically the 
most relevant process for the most degenerate cases, depends only on the strong coupling and 
rrif and does not depend on y x or A h- Finally, the leading contribution to the annihilation of 
tt* pairs into pairs of massive electroweak bosons are controlled by the interactions of t with 
the Higgs boson and therefore the annihilation rates scale as \\. Typically these processes 
are subdominant compared to the annihilations into gluons, however, they can become more 
relevant once A^ > 1 is considered. 
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B. QCD Sommerfeld Effect 


The annihilation of two scalar top partners, which plays a pivotal role in the determination 
of the dark matter relic density, is affected by the non-pertnrbative Sommerfeld effect. As 
a result, the annihilation cross-sections at the time of freeze-out in the different final states 
can significantly differ from the Born-level values given in Table [T] We briefly review here 
the formalism to include the Sommerfeld enhancement in the calculation, following closely 
the approach of Ref. mt 

We consider an annihilation process where the particles in the initial state interact with 
each other via a long-range interaction described by a central potential. For a cross section 
with partial wave expansion a = Xu=o a i y2l ~ 1 ^ the Sommerfeld corrected cross-section is 
given by 

v s = aiv 2l ~ l , (14) 

1=0 

namely, the non-perturbative effects can be calculated separately for each partial wave and 
are encoded in the enhancement factor 

Interactions mediated by a massless particle generate a Coulomb-like potential of the 
form V = a/r, attractive for a < 0, where a is the potential strength. The enhancement 
factor due to such a potential for the s-wave component of the annihilation cross section is 
given by 

_ —na/f3 
0 i _ e W/3 ’ 



while for l > 0 


Si>o = <So X n (1 + 


n= 1 


a 


4/3 2 7i 2 


(16) 


where /3 — v/2 (recall that v is the relative velocity of the incoming particles). 

It follows from this expression that when v«l, S)>o ~ SqV~ 21 , which seems to jeopardize 


the convergence of the partial wave expansion, Eq. (14). However, we will be interested in 
potentials where a ~ 0.1, therefore the terms with high l are suppressed by a factor (0.1 21 /l\ 2 ) 
which ensures the convergence of the series. This implies that, in practice, it suffices to keep 
only the first few terms in the expansion. Then, casting the Born-level cross section as 


a 


- + a 1 v + 6 l , 
v 


(17) 
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where 81 includes the higher l contributions to the cross-section, the Sommerfcld corrected 
cross-section is well approximated by 


a s ~ S 0 — + Sia-L v + S 2 Si , 
v 


(18) 


which we implement in our numerical work. 

The annihilations (it* —>■ hh,\\V 2 ), with V\, V 2 being gauge bosons, can be significantly 
affected due to the potential generated by gluon exchange. Prior to the QCD phase transi¬ 
tion, this potential can be approximately described by a Coulomb-like potential p5| 36| 


V(r) « C— 
r 


a s (li = l/r) 
2 r 


[Cq — Cr — Cr>] 


(19) 


where a s is the strong coupling constant evaluated at the scale (i — l/r, Cq is the quadratic 
Casimir coefficient of the color representation of the final state and Cr and Cri are the 
quadratic Casimir coefficients of the incoming particles. In the case of scalar top partners, 
Cr = Cri = C3 = Cg = 4/3, whereas Cq depends on the annihilation channel. The relevant 
final states are characterized by C\ = 0 for a singlet and Cg = 3 for an octet representation. 
As the sign of C depends on Cq, the potential can be both attractive or repulsive, depending 
on the color representation of the final state. Consequently, annihilations into color singlet 
final states, e.g., hh, receive a universal enhancement, whereas the case is more complicated 
for annihilations into final state which can be in more than one color representation. 

An initial it* state is decomposed as 3 ® 3 = 1 © 8. 1 The QCD potential corresponding 
to these states is 

a, f —4/3 for the 1, 

V = — x { ' (20) 

r y +1/6 for the 8. 

The most relevant final states for cosmology resulting from such an initial state are gg, hh, 
W + W~ and ZZ. For completeness we also include the 'yg and Zg final states. Only the 
gg final state can have both the singlet and octet states, the others are either only color 
singlets or color octets. For the gg final state the Sonmierfeld enhancement factor for each 
partial wave reads mi 


g tt*->gg 



+ 


5 

7 


Si 



( 21 ) 


1 A tt initial state, which decomposes as a 3 ® 6 can be relevant for larger masses than those considered 
here. We include it, and the relevant Sommerfeld enhancement, in our numerics. 
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where the coefficients can be derived from a simple group theory calculation na. Whereas 
for the singlet/octet only final states, 


5S 


ti*^hh.,ZZ,W+W~ 


= s. 


gtt*-tZg,'yg _ g 


a=— a a a 


( 22 ) 


with Si as given in Eqs. (15) and (16). 

The above description of the Sommerfeld effect holds if a definite QCD representation 
can be assigned to the initial state. It has been argued in Refs. [SUEZ! that this assumption 
might not be correct in the thermal bath present at freeze-out: rapid interactions with 
gluons could continuously change the color state of the initial particles and thus prevent 
the formation of a definite QCD color state. In this case, a color averaged initial state 
could be considered instead. Qualitatively, this prescription corresponds to the replacement 


of the individual QCD potentials of different representations in Eq. (19) with an averaged 
potential. In that case all the cross section with the same initial state should be scaled by a 
single Sommerfeld factor independently of the final state. Performing the color average, we 
find the following QCD potential 


T t 1 11 

V ~ X _ 42 


(23) 


We will refer to the above as the color averaged Sommerfeld effect. We explore the effect 
of these two prescriptions in what follows, but primarily use the “color coherent” prescrip¬ 
tion. Regardless of the prescription chosen, the picture remains qualitatively the same. A 
definitive answer to the proper treatment of these thermal effects lies beyond the scope of 
this work. An additional thermal effect arises from the screening of the QCD potential due 
to the gluon’s plasma mass. As was shown in Ref. 021. including the thermal mass of the 
gluon has an imperceptible impact on the relic density. 

Finally, it should be kept in mind that dark matter annihilation and the Sommerfeld effect 
are processes happening at different scales. Consequently, it is not appropriate to evaluate 
the strong coupling a s {/i) which enters the enhancement factors Si at the energy scale of 
the hard annihilation process // m 2 m x . Rather, the Sommerfeld scale fi ~ p, where p is 
the momentum of the annihilating particles, should be used. This could bring into question 
the validity of our calculation if a s (/i) enters the non-perturbative regime. However, this is 
not an issue because the low abundance of such very small momentum particles renders the 
details of the running of a s at these scales unimportant. 
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m x (TeV) 

(a) 



m x {TeV) 
(b) 


FIG. 1: a) The mass difference, Am, required to obtain the correct relic density with (green) and 
without (dashed red) the Sommerfeld effect. The green band spans the mass region between the 
color averaged (dashed) and the color summed (solid) calculation of the Sommerfeld enhancement, 
b) The percent contribution of the various annihilation channels to the relic density, varying Am 
as needed for each m x to satisfy Qh 2 = 0.12. 

C. Numerical Results 


We implemented our simplified model using CalcHEP 3.4 [32] in MicrOmegas 3.3 
and incorporated the appropriate Sommerfeld factors for the different annihilation channels 
relevant for coannihilation of a scalar top partner as described above. We now discuss the 
cosmology of this model. Experimental constraints and future prospects on the presented 


parameters due to collider and other astrophysical observations will be discussed in Sec. IV 


1. MSSM: B + t R 

The MSSM limit of our simplified model is reached by identifying x with a bino like 
neutralino B and the scalar t is specifically a right handed stop t R . The second stop is heavy 
and therefore decoupled. In terms of the MSSM parameters, the demand that the stop mass 
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eigenstates be the gauge eigenstates, explicitly requires a vanishing mixing parameter in 
the stop sector, X t = 0. For consistency with the observed Higgs boson mass, the heavy 
stop would then need to be PeV scale. The couplings in this case approximately reduce to: 


Vx 


MSSM 


= | V2g Y and X h = y 2 . 


In Fig. la we show the mass splitting required between the B and the tn to accommodate 


the observed relic density. The red dashed curve denotes values neglecting the Sommerfeld 
effect. The band of green values is obtained by considering the two prescriptions for the 


Sommerfeld enhancement factor, as detailed in Sec. Ill B 


Fig. lb shows the annihilation channels giving the dominant contributions to the relic 
density as a function of the B mass. We have not displayed (tt* —y Zg/^g), which combine 
to contribute % 10%. The importance of the different channels depends strongly on m x . 
The channel (tt* —> gg) is only dominant when m x is greater than about 700 GeV. If the 
relic density is computed just using this channel, one would obtain a mass difference approx¬ 
imately 10 GeV smaller than the results shown here HUE!- Again the width of the band 
corresponds to the different treatment of the Sommerfeld prescriptions. For the remainder 
of the paper, we will utilize the color coherent prescription, bearing this uncertainty in mind. 


2. Simplified Model 

In our simplified effective stop model, the couplings y x and A h are free parameters. The 
quartic coupling A h only impacts the coannihilation channels which are related to the Higgs 
boson, either in the intermediate or the final state. 

As a means of analyzing this multidimensional parameter space we will first fix Xh to a 
given value. Then, we find the minimum r = r min that can yield the correct abundance, 
found by saturating the necessary annihilation cross section by the gg final state alone. For 
all values of r > r m i n there exists a y x such that the experimentally observed relic density 
can be found. In Fig. [2a] and [2b] we show the contributions to the relic density from different 
initial states in the m x - r plane for Xh = y 2 ~ 1 (left) and A h ~ « 4 (right) respectively, 

such that the relic density is fixed to be filh 2 = 0.12 by varying y x . The white portion of the 
plots denotes the region where the experimentally observed relic density cannot be obtained 
thermally, even with y x = 0. The red, bright green and blue regions denote where the 
dominant initial states (i.e. contribute more than 50%) are tt, \t and XX respectively. In 
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m x ( GeV ) nty ( GeV ) 

(a) (b) 


FIG. 2: Regions where the relic density is dominated by different initial states, a) A/,, = yf ~ 1. 
The black line correspond to the minimal MSSM B + Ir model, b) A/,, = 4 yf ~ 4. T/ie black line 
corresponds to y^ ISSM = |\/2gy. 


the light green shaded region several of the three initial states are of similar strength and no 
dominant channel can be identified. The black line in both plots shows where y x is equal to 
the MSSM value of 2y/2gy/3. Comparing the left and right panels we can see the impact 
of changing the value of the quartic coupling A h on cosmology. Defining coannhilation as 
the region where the yy initial state is not the dominant contribution, we can see that, 
apart from the very low mass region, the value of r for the cross-over from coannihilation to 
annihilation (i.e. the transition to the blue region) for a given m x , is relatively insensitive 
to A h- Large A/,, does effect the coannihilation region: For A^ ~ 4 the red region shifts 
approximately 3% upwards in r across the m x range considered with respect to the A^ ~ 1 
case. At very low masses, there is a change in this behavior. This can be understood by 
comparing the A h dependence of the three channels, it* —* hh, WW, ZZ: At lower masses 
of t the scaling of the hh channel transitions from \ 2 h to A^ (cf. Table [ij. In any event, as 
we will see, a light scalar top partner with such a large quartic coupling to Higgs bosons 
would cause significant deviation in gluon fusion, which is most certainly excluded by Higgs 
coupling measurements [38]. 
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n h 2 = 0.12 



FIG. 3: y x as a function of r = m^/m x for various choices of the dark matter mass. The shaded 
regions corresponding to each mass are bounded by Xh = y\ ! ~ 1 (solid) and Xh = 4 y 2 ~ 4 (dotted). 
The A h = 2 y( ~ 2 value is denoted by the dashed lines in each shaded region. The horizontal black 
dot-dashed line denotes the value of the MSSM coupling, y x SSM . 


Fig. [3] shows values of y x required as a function of r for different m x to obtain flh 2 = 0.12. 
The width of the colored bands captures the effect of varying A^, with the solid (dotted) lines 
corresponding to A^ ~ 1 (4). The A^ ~ 2 value is shown by the dashed line in the middle 
of the shaded regions. The horizontal black dot-dashed line denotes y x = y^ rsSM = | \f2gy, 
corresponding to the MSSM value. The values of r where the y(\ ISSAI line intersects the 
A/j ~ 1 lines for each m x corresponds to the mass splittings shown as the solid green curve 


in Fig. la The narrowing of the colored bands for increasing r shows the transition from 
coannihilation into the self-annihilation region and confirms that the value of Xh is only 
relevant for cosmology when there is a significant contribution to the relic density from 
coannihilation. Further, as also seen from the previous figures, even extreme values for A h 
only result in a shift of at most 3% for the needed r (for fixed m x ) throughout the mass 
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r = m~ t /m x 

FIG. 4: Percentage contributions from relevant annihilation channels as a function of r = mi/m x , 
resulting in Flh? = 0.12 for m x = 500 GeV . The shaded bands denote the variation in A h from 1 
to 2. The corresponding y x is as shown in Fig. [3| 

range we consider. In contrast, values of y x comfortably below perturbativity constraints 
can allow for large mass splittings, easily allowing for r > 1.2 even for the heaviest m x 
considered. 

In Fig. [4] we show the relative percentage contributions from different channels to the 
total relic density, given Vth 2 = 0.12, as a function of r for m x = 500 GeV. The bands 
correspond to the variation of A^ between 1 (solid lines) and 2 (dashed lines). The values of 
y x correspond to those shown in Fig. [3j As A h increases, the smallest value of r consistent 
with the proper relic density increases. Note that since the contribution due to the ( tt* —>- gg ) 
channel depends only on the mass splitting for a given m x , it doesn’t get perturbed by the 
change in A*,. Therefore, the increase in (tt* — y hh/VV) is compensated by a decrease in all 
the other channels shown. 
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IV. EXPERIMENTAL CONSTRAINTS AND FUTURE PROSPECTS 


Having determined the regions of parameter space in which thermal freeze-out can account 
for the observed relic density, we now turn to other searches for dark matter. We will discuss 
briefly the relevant characteristics of the different experiments before analyzing the impact 
of these observational efforts on the parameter space of thermal dark matter. 

A. LHC 

Dark matter that primarily couples to top quarks suffers from a severely suppressed direct 
pair production cross section. Therefore, the most promising way to probe this type of model 
is via the production of the mediator t. Scalar top partners are charged under QCD, so pairs 
of the mediator can be produced copiously at the LHC, with the rate depending only on the 
mass mp The signature of t depends on the available decay channels which are primarily 
governed by the mass difference between the dark matter and the t. As long as m^—m x > mt , 
the decay (t —>• yf) will be dominant, whereas in the range m t > rrif — m x > rri\y + mb, the 
three-body decay, (t —> yW6), becomes relevant. Therefore, the LHC limits on our model 
are similar to those on models for direct stop production with branching ratios of 100% into 
these channels and the limits derived by the ATLAS pITJl !TT| and CMS jT?i 155] collaborations 
apply without any modifications. However, when my/ + mj > — m x , the region in which 

we are predominantly interested, the flavor-violating process, (t —> yc), as well as the four 
body decay, (t —> y bff), could in principle contribute significantly. In our model, absent 
additional flavor violation, the four body channel would dominate. However, various UV 
completions to our simplified model could allow the relative contribution of these channels 
to be a free parameter without changing the cosmology. Within the MSSM, the branching 
ratio into both the decay modes can be significant without violating flavor constraints, see 
e.g., Refs. [35, 136] . 

Both ATLAS and CMS have analyses optimized for each of these channels, under the 
assumption of 100% branching ratio into either decay mode [57,(58]. At small mass splittings 
the monojet search is interpreted for both the channels, and found to have similar sensitivity. 
This is unsurprising; at very small mass splittings, one would expect the final state particles 
to go undetected, for either of the two decay channels. As the mass splitting increases, 
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particles in the four body decay might run afoul of lepton or jet vetoes in the present 
analysis. A dedicated search with charm-tagging is relevant for the 2 body decay, and has 
sensitivity primarily when the mass splitting is larger than about 20 GeV B7j. Also at 
larger mass splitting (but still less than mw) if the 4-body decay dominates, there is some 
sensitivity to final states using soft leptons B3- In this moderate mass splitting regime, it 
has been emphasized that the exclusion limits are in fact very sensitive to this branching 
ratio and can be weakened considerably when both channels have competing branching 
ratios [49]. In our analysis, we impose limits obtained from ATLAS, similar results would 
be obtained if the CMS limits were used instead. We will display both a monojet region 
(which we expect to be insensitive to the branching ratio), and a region that explicitly relies 
on charm tagging. 

In addition, with the given mass splittings between the bino and the stop, even though 
LHC14 can only probe up to ~ 500 GeV in stop masses, the entire mass range consistent 
with the observed relic density can be comfortably probed by a 100 TeV collider [50j. These 
projections are for monojet searches, without any assumption about the decay mode, and 
so should be quite robust. 2 

There can also be a significant impact on gluon fusion for the observed 125 GeV Higgs 
boson due to the presence of t. The gluon fusion amplitude can be simply written down 
using the low energy theorem [5TH53] : 


A , ~ ^sm , 

•^hgg — S^hgg ' 


A hV 2 

2 m 


2 > 


(24) 


where we have used the normalization that the SM contribution due to the top loop is 
*4.^=4. 3 The above is an excellent approximation in the limit that the relevant t mass 
is sufficiently heavier than the Higgs boson mass. The best fit signal strength reported by 
CMS for the 125 GeV Higgs boson from gluon fusion is Hh gg = 0.85^o;)g [3H] • Since scalar 
top partners without mixing can only give an enhancement, it is reasonable to impose the 


2 We reemphasize, that when the mass splitting approaches mw, the sensitivity of the monojet search for 
the four-body decay interpretation will be reduced if the lepton and jet vetoes employed by the ATLAS 
and CMS collaborations for the current monojet searches are used. 

3 A similar contribution is induced in the diphoton decay width of the Higgs due to the presence of 

t. However, the SM contribution due to the W and t loops has opposite sign and is much larger in 
magnitude, = —13 [511156] . Additionally, since the rate into 77 is proportional to the dominant 

production mode of the Higgs, gluon fusion, the total impact of a scalar top partner on the diphoton 
decay rate is further diluted. 
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requirement that the contribution of t to gluon fusion does not exceed the SM value by more 
than 20%. This implies a constraint on mi as a function of A h, which can be rewritten as a 
constraint on r as a function of m x and A h- 


> 


277U 


A h 


2 {y/f'hgg l) 


1/2 


1-14V%i —, 


77U 


(25) 


assuming fihgg is bounded to be less than 1.2. Clearly, as experimental precision and the¬ 
oretical control increase, the above constraint will become stronger and, as we will see, is 
already quite restrictive for the case of large A h- 


B. Indirect Detection 


The annihilations (xx ~t tt) are expected to occur today in regions with an overdensity 
of dark matter particles, thus leading to potentially observable signals in indirect detection 
experiments. In these regions, dark matter particles typically have very low velocities, 
therefore, the self-annihilation cross section is almost entirely s-wave, in contrast to the time 
of freeze-out, where both the s- and p- wave contribution were relevant. Using Eq. (§. in 
the absence of coannihilation, the annihilation cross section today for a thermally produced 
dark matter particle can be approximated by 

av ~ ( -~ 1/2 3 x 10- 26 cm 3 /s , (26) 

\a + 3b/x F J V25/V80/ ' K ’ 

where a and b are given in Eqs. © and 0- respectively, and x F ~ 25. In particular, 
we find that crv ranges between (2 — 0.1) x 10 -26 cm 3 /s for m x = 250 GeV — 2 TeV. These 
values of the cross section lie below the present upper limits from H.E. s.s. m and Fermi- 
LAT [58], however, depending on the dark matter density profile in the galactic center, the 
prospected Cherenkov Telescope Array (CTA), might possess the sensitivity necessary to 
probe crv > 1 x 10 _26 cm 3 /s in this mass region [59]. In the coannihilation region, on the 
other hand, the cross section for (xx ~► tt) is highly suppressed, therefore the detection of 
annihilation signals will be very challenging in this regime. 

Recently, after a reevaluation of the background uncertainties and a new derivation of 
the gamma ray spectrum [60], {xX ~► tt) has also been considered as a possible explanation 
of the Galactic Center gamma-ray excess [611162]. Interestingly, it was found that for m 


X ~ 


200 GeV, annihilation of thermal bino dark matter into top pairs might account for the 
observed excess, even though the p-value is rather low 
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(a) 


(b) 



FIG. 5: Representative examples of the triangle ((a) and (b)) and (c) box diagrams which 

contribute to the dark matter nucleon coupling. 


C. Direct Detection 


At first sight the prospects for the observation of dark matter coupling to top quarks 
do not seem to be particularly encouraging as the absence of top quarks in the nucleus 
prevents tree-level interactions. However, we find that loop diagrams can change this picture 
considerably and can induce a spin-independent dark matter nucleus scattering cross section 
within the reach of upcoming experiments. 

The first process which can generate a sizable direct detection cross section is due to a 
loop induced coupling of the dark matter with the Higgs boson. The effective xx) 1 coupling 
is generated by triangle diagrams with scalar top partners and top quarks in the loop, 
see Figs. and it for diagrams. This effective coupling has been calculated in the case 
of supersymmetry |63| and the results can be applied to our model with the appropriate 
replacements. We re-derived the dark matter Higgs boson coupling using the low energy 
Higgs theorem [FI] and find good agreement between our calculation and the result of 
Ref. [63], Loop induced interactions with the Z, which are known to be important for 
Dirac dark matter interacting with tops [63], are not expected to be relevant here as vector 
interactions vanish for Majorana fermions. 

The second relevant loop effect is due to the fact that the dark matter can couple to the 
gluon content in the nucleus via a box diagram with scalar mediators and top quarks in the 
loop (e.g. Fig. [ 5 J 3 ) |65] . 

Another effect that cannot be neglected when A h is large, is due to the loop induced 
coupling of the Higgs boson to gluons via t. However, this will be a two-loop effect in the 
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direct dark matter detection cross-section. Furthermore, this is precisely the loop which leads 
to modifications of the gg —y h rate, and hence the region where this would have a relevant 
effect on the direct dark matter detection cross-section would lead to an unacceptably large 
deviation in the production of the Higgs boson at the LHC. Therefore, we do not consider this 
effect in our analytics below. However, it is always included in our full numerical calculations 
which were computed using the built-in analytical MSSM formulae in MicrOmegas 3.3, with 
appropriately rescaled couplings, including the effective yy/i coupling computed according 

to Ref. [63], 

The direct detection cross-section per nucleon is given by: 


Cor — 


Arri^rri^ A 2 


(27) 


s/ n(m x + m„) 2 ’ 

where n refers to either the neutron or the proton and A is the amplitude. A can be written 
as follows: 


A = A h + A », 


(28) 


and we will further decompose the contributions in the amplitude due to the exchange of 
the Higgs boson in a part which depends on A h (where the t couples to the Higgs boson, 
Fig. |5^) and independent of A h (where the t couples to the Higgs boson, Fig. [5 ]d): 


A h 


W+Ailvf” 

2 vmi ^ q ’ 

n n 


frq, Q = { u , d , s} 

§ffr g , Q = {c,b,t} 


(29) 


where = (1 - fr q ) and f Tq = {0.0153,0.0191,0.0447}, f? q = 

{0.0110, 0.0273, 0.0447} are the default values used in MicrOmegas 3.3, leading to E q fq — 
0.28. 


As the full analytic expressions for both the effects described above are lengthy and 
cumbersome we do not repeat them here and refer the reader to Refs. (63j [65J. However, 
certain expansions can be made in the parameter regions of interest which gain us insight 
in the behavior of the various contributions. 


From our numerical results in Sec. Ill C we know that the mass splitting required in 
the MSSM to obtain the experimentally consistent relic density is less than ~ 45 GeV, 
for m x < 2 TeV. When we allow y x and A h to be free parameters, the mass splitting can 
be significantly larger. In fact, for very heavy dark matter masses with appreciable mass 
splitting, we will see that the direct detection cross-section can be enhanced, partially due 
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to the fact that a large y x is needed to obtain a consistent relic density. Therefore, we will 
present analytical expressions for two regimes: Am smaller or larger than the top mass. 

When the mass splitting between m x and is smaller than we can expand in the 
small parameters S = {m t — mi + m x )/m x and £ = (m t + mi — m x )/m x : 


A ^ 
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3 y x m 
47t 2 vm 
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(31) 

(32) 


where the default value for G l a = 0.053 from MicrOraegas 3.3. Note that when the mass 
splitting, rtii — m x , is much much smaller than m*, then 6 ~ £, however, when the mass 
splitting is close to the top mass, then 5 « £. We have checked the above expressions 
reproduce the full numerical results within 30% for all three amplitudes when m x > 500 GeV, 
leading to an estimation for the total cross-section which is accurate to 50%. For smaller m x , 
the above approximation gives results within a factor of 2 from the full numerical calculation 
as long as r is not too small (6 small). 

Typically one expects the Higgs boson exchange to dominate the dark matter nucleon 


coupling, however, as can be seen from Eqs. (30)-(32), all the different amplitudes scale 
with approximately l/m x in this regime. In particular, this means that the dark matter 
coupling to gluons can become comparable to its coupling to the Higgs boson in this small 
mass splitting regime, especially for dark matter mass, m x , close to m*. Unfortunately, the 
contributions from the triangle and the box diagrams interfere destructively such that direct 
detection experiments have reduced sensitivity to dark matter with m x close to rn t . 

When the mass splitting is larger than the top mass, one can take the limit that the mass 
of the scalar top partner is much larger than both m x and mt. In this case, the various 
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contributions to the amplitude are substantially simplified: 
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(34) 

(35) 


Using the above expansions, the Higgs exchange amplitude A h , defined in Eq. (29), goes as 


1 /{rrihirii) 2 whereas A 9 is proportional to l/m~. Consequently, the contribution from the 
triangle diagrams will always dominate the cross section for large mp Comparing the above 
amplitudes with the full numerical results, the contribution from t, A~, is within 10% of the 
full numerical amplitude, even for small r ~ 1.2 and m x > 500 GeV. The top contribution, 
A%, is accurate to about 20% across the region of interest. The contribution from the box 
diagrams with the gluon A 9 is only accurate to approximately an order of magnitude, but 
the gluon contribution is negligible in this regime. Even ignoring it completely, the direct 
detection cross-section is within a factor of 2 of the full numerical calculation. 

D. Numerical Results 


In the following we discuss the impact of various experimental probes described above 
on the allowed parameter space of thermal dark matter. We begin with the description of 
the effect on the more general simplified model. We also draw conclusions for our minimal 
MSSM scenario. 

As discussed in the previous subsection, in the simplified effective stop model, the loop 
induced couplings between the dark matter and the SM can have a significant impact on the 


direct detection cross-section. In addition, as pointed out in Sec. IV B in the blue shaded 


region shown in Figs. 2a and 2b, annihilations could give rise to detectable signals at future 
indirect detection experiments [59]. The measurement of the gluon fusion cross section at 


the LHC is sensitive to the value A/, and excludes a large portion of the parameter space 
under consideration for the largest value we consider. We stress again that direct searches 
at the LHC do not depend on either of the couplings, but are only sensitive to the mass 
splitting between the two states and the overall mass scale. 
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FIG. 6: The blue dashed contours label direct detection cross-section in units of ICG 47 cm 2 . Green 
contour lines show the value of y x required to saturate the relic density. The green contour labeled 
“MSSM” corresponds to y x = | y/2gy, and lies below the irreducible neutrino background, denoted 
as “Neutrino n w- The bright green shaded region corresponds to the region which will be probed 
by LZ JE|/. The region not shaded green is where it is not possible to saturate the relic density 
constraint thermally in this scenario. Also displayed are lines showing where the mass splitting 
Am = ruf — m x = m\y or mt- Red (dark: monojet, bright: charm-tagging) and orange (monojet) 
regions denote current 139 1 \2fl\) and projected f exclusion bounds from the LHC. The region to 
the left of the black line denoted as “Gluon Fusion” would give rise to more than 20% enhancement 


in gluon fusion compared to the SM expectation, Eq. (25). 


Figs. [6] and [7] summarize our results for the A/, - 1. 2 and 4 cases respectively. Green 
contour lines show the value of y x required to saturate the relic density. The green contour 
labeled “MSSM” in Fig. [6] corresponds to the minimal supersymmetric model for our set up, 

given by y x = \ > \[ 7 2g Y (also denoted in Fig. [t] as the contour labeled 


discussed in Sec. 


Ill Cl 


.MSSM \ 


y x """). The region not shaded green is where it is not possible to saturate the relic den- 


24 













2000 


200 300 500 700 1000 



(a) 


200 300 500 700 1000 2000 



200 300 500 700 1000 2000 

m x (GeV) 


(b) 


FIG. 7: Same as Fig. but showing results for (a) Xh = 2~ 2 and (b) A h = ~ 

shaded region in (b) corresponds to the region that will be probed by XEN0N1T ]66f . 
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sity constraint thermally in this scenario. The blue dashed contours label direct detection 
cross-section in units of 10~ 4 ' cm 2 . As a comparison we indicate the so called neutrino floor, 
labeled as “Neutrino”, where the rate of dark matter interactions in a detector equals the rate 
of recoil events induced by coherent scattering of neutrinos (primarily atmospheric neutrinos 
for the masses we consider). As these neutrino interactions constitute an irreducible back¬ 
ground, direct detection experiments cannot be expected to probe appreciably smaller cross 
sections. This irreducible background interpolates between ~4x 1CT 49 — 4 x 10 -48 cm 2 for 
the mass range considered [SB]. Except for the X h ~ 4 case, the direct detection cross-section 
for y^ssM lies mostly below the neutrino background floor. The bright green shaded region 
corresponds to the region which will be probed by LZ [66]. XENON1T will have sensitivity 
in the A^ ~ 4 case, which is depicted by the shaded blue region [66]. 

Also shown are lines indicating where the mass splitting Am = mi — rn x = m\y or m t . 
Darker red and red regions denote current LHC exclusion regions due to monojet and charm¬ 
tagging respectively PHUT]. The orange region denotes the projected LHC14 monojet ex¬ 
clusions [50]. Recall that the exclusions due to the monojet searches are expected to be 
model-independent; however, the exclusion limits obtained via charm-tagging could be sig- 
nihcantly weakened as discussed in Sec. IV A The region to the left of the black line denoted 
as “Gluon Fusion” would give rise to more than 20% enhancement in gluon fusion compared 
to the SM expectation, Eq. (25). Outside of the region plotted in Figs. [6] and [7] are additional 
LHC constraints for Am > rn t in the low m x , large r region. Specihcally, for Am ~ m t , 
even though current bounds only extend to m x ~ 275 GeV [401 [UJ ITT] , LHC14 is expected 
to probe m x ~ 500 GeV m- These searches will not probe the coannihilation region, but 
will be complemented by possible signatures in future indirect detection experiments [ 59] . 

We now discuss the direct detection contours in some detail. We recall that the direct 
detection cross-section is always proportional to y x , which we Sx by requiring a consistent 
relic density at every point. This requirement implies that for every m x there is a max¬ 
imum occurring at some r and a si is increasing with m x . To see this we Srst note that 
coannihilations begin to suppress the required y x severely for r < 1.1 (c.f. Fig. [3]), and so 
the smallest direct detection cross-sections are found in the most degenerate region. The 
maximum occurs approximately when the coannihilation processes become irrelevant for 
setting the relic density, and the relevant process is (xx tt) (corresponding to when the 
contribution of this channel is > 70%). This occurs for r ~ 1.15 — 1.3 across the mass range 
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we consider. For larger r, y 4 can be determined by examining the partial wave expansion 
for the (xx tt) process. In particular, y 4 is inversely proportional to (a + 3 b/xp), with 
values set as in Eqs. (10) and 




3 m; 
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3 (r 4 + 1) 




327T (r 2 + 1 - nril/mty m 4 167r (r 2 + 1) m 2 x F 


(36) 


One expects that in the region where the s-wave contribution is dominant, p 4 scales as m 4 
and as m 2 in the region where the p-wave contribution is the most relevant. Due to the 
different r dependence of the s and p wave contributions, the annihilation cross-section is 
not dominantly p-wave until almost m x ~ 2 TeV. 

Now turning to the scaling of asi with m x and r, we note that in the small mass splitting 


regime, Eqs. (30)-(32), 
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Instead in the large mass splitting regime, Eqs. (33)-(35), 
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Comparing Eqs. (36)-(38), we see clearly that for a fixed r, if the cross-section setting the 
relic density is predominantly s-wave, then the direct detection cross-section increases with 
m x : a si oc m 2 . When the relic density is instead set by the p- wave contribution, the 
dominant scaling of the annihilation cross-section and the direct detection cross-section are 
the same, and we expect very little sensitivity to increasing either m x or y x . Consequently, 
the direct detection cross section of thermal dark matter increases with m x in the mass 
range of interest. 

The stronger variation with r of a si for increasing m x can also be understood by com¬ 
paring the r dependance in Eq. (37), using the first term for p 4 in Eq. (36) (small mass 
splitting and s-wave dominated annihilation cross-section) with the r dependance obtained 


in Eq. (38) using the second term for p 4 in Eq. (36) (large mass splitting and p-wave domi¬ 
nated annihilation cross-section). 

We also comment briefly on the scaling with A/p Since the leading dependance of both A\ 
and A~ with m x and r is approximately the same and A g is only relevant in small regions 
of parameter space, the direct detection cross-section approximately behaves as (p 2 + A^) 2 . 
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It is interesting to note the complementarity in the reach of the LHC and direct detection 
experiments. As mentioned previously, the direct search bounds from the LHC do not 
depend on the exact value of the couplings, but only on the mass splitting and the mass 
scale. The LHC14 should be able to probe masses up to 500 GeV in the Am < mw region, 
a 100 TeV collider would be able to comfortably probe masses up to ~ 2 TeV [50] . On the 
other hand, both the enhancement in gluon fusion and the direct detection cross-section are 
impacted by the mass splitting and the couplings. For ~ 1, the current direct search 
limits from LHC and bounds from gluon fusion are approximately comparable. The direct 
detection cross-section is very suppressed in the region where the LHC bounds are expected 
to be strongest. On the other hand, for the larger m x , r, region, where LHC searches will 
have no sensitivity, the large y x required to obtain an experimentally consistent relic density 
enhances the direct detection rate, which is largest in this region. As can be seen in Fig. [6j 
LZ is expected to probe m x > 600 GeV and Am > mw Once A h is increased to 2, the 
constraint due to the enhancement of gluon fusion becomes stronger, however, the 14 TeV 
LHC is still expected to have stronger sensitivity. The direct detection cross-section also 


increases, and we can see from Fig. 7a that now there is an overlapping region between 
m x = 300 and 500 GeV where both LZ and LHC14 will be sensitive. When A h is pushed 
to an even larger value of 4, the requirement of not having more than a 20% enhancement 
in gluon fusion constrains a large region of parameter space for m x < 500 GeV, somewhat 
stronger than the direct search sensitivity expected from LHC 14. In addition, as can be 


seen in Fig. 7b the direct detection cross-section is enhanced significantly such that certain 
regions will be accessible in the near future to XENON1T. 

Turning now to Fig. [8j we have plotted the maximum achievable dark matter nucleon 
cross section in our model for three different values, A h ~ 1, 2, and 4, as a function of m x . 
For each value of m x , we scan values of rrq and y x consistent with the relic density and 
End the largest cross section. We reiterate that these maximal cross sections are typically 
found where the process (xx tt) starts to dominate the cosmology (coannihilations are 
unimportant). Note that r is relatively constant for the maximal a si values, varying between 
r ~ 1.15 — 1.3, as can be seen from inspecting the direct detection contours in Figs. [6 7b For 
the lower dark matter masses, m x <TeV, as stated earlier, the direct detection cross-section 
increases as However, as the annihilation cross-section transitions to be dominantly 
p-wave, the direct detection amplitude and early universe annihilation both approximately 
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FIG. 8: The maximum value of the dark matter-nucleon cross section for different values of A/,, for 
points consistent with a thermal relic density. Also shown are the expected sensitivities of both the 
XEN0N1T and LZ experiments. We have imposed a maximum value of y x < 3. The exact value 
of this bound is largely irrelevant for determining the maximum cross section, but can effect the 
maximum cross section, particularly at small masses. Moreover, it does determine the maximum 
m x for which a relic density may be successfully obtained. For details, see the text. 

scale as ~ y^/rri^. Thus, once the relic density is fixed, the direct detection cross section 
does not change dramatically even as y x and m x are increased. 

Nevertheless, we do cut off the scan at values of y x < 3 = y™ ax which corresponds to a 
dark matter mass of 3 TeV. For most masses, the value of y™ ax does not affect the maximum 
direct detection cross section - it is determined by cosmology alone. However, at the lowest 
masses (i$ 250 GeV) it would be possible to raise the direct detection cross section very 
modestly (say 10%) by allowing a larger y™ x ~ 5. 

An important take home message is that in WIMP models where the dominant inter¬ 
actions are with the top partners, cosmological considerations indicate a cross section well 
below current direct detection bounds. Indeed, even XENON1T will have difficulty probing 
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much of the parameter space. LZ, however, will be much more effective. 

While we have been considering this simplified model independent of the MSSM, it is 
worthwhile to ask what the above results imply for a MSSM spectrum that mimics the one 
we consider here. In addition to the loop mediated contribution, important contributions 
to direct detection may arise from tree-level Higgs boson exchange which are generated by 
non-zero Higgsino-bino mixing. Assuming that the non-Standard Higgs boson contribution 


decouples, the Higgs boson mediated direct detection of Eq. (29) is modified to 

+ A\ + A tree ) 
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where n is the supersymmeteric Higgs mass parameter, tan (3 = v u /v,j is the ratio of the two 
MSSM Higgs held vacuum expectation values and v 2 + v 2 , — v 2 . 

In the limit where only A tree contributes, the relevant cross section may be approximated 
as: 
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.500 GeV, 

Looking back to Fig. [6j it is easily possible - even for relatively large g, - that this contribution 
can dominate the loop mediated one, which for the MSSM case lies at the 10 -48 cm 2 level 
(Note the cosmology would be relatively insensitive to the presence of a small Higgsino 
admixture.). So in the case of a bino-f# MSSM-like case ( y x ~ 2y/2gy/3, A h ~ 1) an 
observed signal at LZ could be a hint of new dynamics (e.g. Higgsinos) at near the TeV 
scale. Note a more general MSSM with stop mixing could effectively allow larger stop-Higgs 
couplings (due to large A-terms), which can modify the direct detection and cosmology in 
important ways. This is the subject of upcoming work. 


V. CONCLUSION 

We have investigated a simplified model of a dark matter candidate to gain insight into 
models where interactions with a top partner are important. This applies, for example, to 
stop coannihilation in the MSSM. 

We showed that for the supersymmetric B coannihilating with the tn, in the low mass 
region, m x < 500 GeV, channels other than the final state gg can be important and that 
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ignoring these channels can lead to a shift of approximately 10 GeV in the expected mass 
splitting between m x and mf . The LHC has the potential to cover the coannihilation sce¬ 
nario np to m x ~ 500 GeV, and a 100 TeV collider can cover the entire range of dark matter 
masses consistent with the relic density, m x < 2 TeV. Unfortunately, the expected rate of 
dark matter interactions in direct detection experiments remains too low for a discovery 
of this minimal supersymmetric scenario, even when loop mediated processes are consid¬ 
ered. Confirming the dark matter interpretation of a possible collider signal by additional 
astrophysical observations is going to be very challenging. Conversely, a signal in any di¬ 
rect detection experiment can only be accommodated within this simple MSSM scenario if 
additional new physics is at the TeV scale. 

Allowing for a free y x , one can determine the required coupling to obtain a consistent relic- 
density for any given mass splitting between y and t as a function of the dark matter mass, 
m x . We note that A h is much less powerful than y x for cosmology and only relevant in the 
small r = m^/m x , region. However, can have a significant impact on gluon fusion, which 
is the main production mode for the 125 GeV Higgs boson. Collider searches are insensitive 
to the couplings y x and A h, but depend on the mass splittings, therefore a large region of the 
parameter space under consideration will be probed by LHC 14 and, more comprehensively, 
at a 100 TeV collider. 

We also computed the loop induced coupling of a pair of x to the Higgs boson, which 
is particularly relevant for the direct detection cross-section. The direct detection cross- 
section scales with y* and depends on m x , r and A/,. While current bounds from LLIX are 
not able to constrain this scenario, we find that the sensitivity of near future direct detection 
experiments, namely LZ and XENON1T, will allow for the testing of this scenario. 

In conclusion, we showed that our Simplified Model has a rich and interesting phenomenol¬ 
ogy. Current experimental limits leave much of the parameter space untested. However, it 
is interesting to note that a combination of future collider direct searches and indirect and 
direct dark matter detection experiments will comprehensively probe the parameters of this 
model. 
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